Stochastic Model for a Vortex Depinning in Random Media 
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We present a self-organized stochastic model for the dynamics of a single flux line in random 
media. The dynamics for the flux line in the longitudinal and the transversal direction to an averaged 
moving direction are coupled to each other. The roughness exponents of the flux line are measured 
for each direction, which are close to ay ~ 0.63 for the longitudinal and a_i_ « 0.5 for the transversal 
direction, respectively. The dynamic exponents are obtained as 2 ~ 1 for both directions. We 
discuss the classification of universality for the stochastic model. 
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In the past few years, there has been an explosion of 
studies in the field of dynamics of fluctuating interfaces 
due to theoretical interests in the classification of univer- 
sality for stochastic models and also due to applications 
to various physical phenomena such as crystal growth, 
vapor deposition, electroplating, biological growth, etc. 
A number of discrete models and continuum equations 
for interface dynamics have been introduced and studied 
[1-3]. An interesting feature of nonequilibrium interface 
dynamics is the nontrivial dynamic scaling behavior ^] 
of the interface fluctuation width, i.e., 



WiL,t) ^ Y.^h{x,t) - m)Y'^ ~ L^it/L'^), 



(1) 



where h{x,t) is the height of site x on substrate at time 
t. h, L, and d' denote the mean height, system size, and 
substrate dimension, respectively. The angular brack- 
ets stand for statistical average. The scaling function 
behaves as f{x) — > const for x 1, and f{x) ~ x'^ 
for X <^ 1 with z = a/(3. The exponents a, f3 and z 
are called the roughness, the growth, and the dynamic 
exponents, respectively. 

Recently the problem of the pinning-depinning (PD) 
transition of interfaces in random media has also at- 
tracted interests in association with the problem of dy- 
namics of fluctuating interfaces. Examples include the 
dynamics of domain boundaries of random Ising spin 
systems after being quenched below the critical temper- 
ature 1^, wetting immiscible displacement of one fluid 
by another in a porous medium ||^, pinning flux lines in 
type-II superconductors 0| , fluid imbibition in paper , 
etc. In the problem of the PD transition, interface is 
pinned when external driving force F is weaker than pin- 
ning strength induced by random media, while it moves 
with a certain velocity v when the force F is greater than 
the pinning strength. Thus there exists a threshold of ex- 
ternal applying force Fc across which the PD transition 
occurs. The role of the order parameter is played by the 
mean velocity, v ~ {J2x^^i^''^)/'^^)/-^'^ ■ Accordingly, 



the velocity is zero for F < Fc, and increases for F > F^ 
as 

v-iF~Fcf, (2) 
where the exponent 9 is called the velocity exponent. 

The continuum equation for the dynamics of interfaces 
in random media may be written simply as j^] 



dh{x, t) 
dt 



i'W^h + F + T]{x,h), 



(3) 



where h{x, t) is the height of the interface at position x at 
time t. The first term on the right-hand side is from the 
smoothening effect of surface tension, the second term 
the uniform driving force, and the third a random force 
with short range correlations, satisfying {ri{x, h)) — 
and {rj{x, h)r]{x' , h')) = 2D6{x — x')6{h — h') with noise 
strength D. The above equation, called the quenched 
Edwards- Wilkinson (QEW) equation, would be relevant 
to the dynamics of the domain wall in random magnetic 
systems. More generally, recently a new continuum equa- 
tion was introduced |lO| , which includes a nonlinear term 
■|(V/i)^ induced from the anisotropic property of the pin- 
ning strength. Thus the equation is replaced by 



dh[x, t) 
dt 



^iyV^h+^{Vhf + F + T^{x,h), (4) 



which is called the quenched Kardar Parisi Zhang 
(QKPZ) equation. The QKPZ equation leads to a differ- 
ent universality class from the QEW equation. Recently 
several stochastic models in the QKPZ universality class 
have been introduced [9,11]. From the models, it has 
been naturally concluded that the surface at the thresh- 
old of the PD transition Fc can be described by the 
directed percolation (DP) cluster spanned perpendicu- 
larly to the surface growth direction in 1-1-1 dimensions. 
The roughness exponent a of the interface is given as 
the ratio of the correlation length exponents and 
of the DP cluster in the transversal and the longitudinal 
direction that is a = v^/u^ « 0.63. 

The dynamics of a single flux line in type-II supercon- 
ductor with random impurities can also be understood 



in a similar manner used in the dynamics of surface 
growth. The main difference between them hes in that 
the flux line is an one-dimensional chain embedded in 
three dimensions rather than in two dimensions. Thus 
the roughness of the flux line is quantified in two dif- 
ferent directions, the longitudinal and the transversal 
to the direction of its averaged velocity. Recently the 
continuum equations for the flux line dynamics in each 
directions were derived by Ertas and Kardar 0, which 
are coupled each other and look very complicated in 
general. They obtained the roughness and the dynamic 
exponents for various cases of the coupled equations; 
however, there still remain several cases that the rough- 
ness and dynamic exponents are not determined. In this 
paper, we will introduce a simple self-organized stochas- 
tic model, which we believe is relevant to the dynamics 
of the flux line. The numerical results we obtain from 
the stochastic model will compensate for the roughness 
and the dynamic exponents for the unknown cases. 

The stochastic model we introduce in this paper is 
defined as follows. First, we consider a body centered 
cubic (bcc) lattice, in which an elastic string runs along 
the x-direction as shown in Fig. 1. The discrete version 
of the elastic string is composed of L-massless beads 
(black dots) which locate at nearest neighboring sites of 
one another and they are connected with strings. The 
configuration looking a zig-zag as shown in Fig.l is re- 
garded as the initial flat configuration. The total length 
of the string is equal to \/2iaL/2, where a is a unit lattice 
constant of the bcc lattice and L is the total number 
of beads in the system. The elastic string does not run 
backwardly to the x-direction, so that (y, z) positions of 
the bead for each x are specified by single values. Each 
bead can be updated only in either positive y or positive 
2:-directions according to the following rule. First, two 
random numbers are assigned on each bead, one of which 
is for positive y-direction and the other is for positive 
z-direction. The random numbers for the ^/-direction are 
uniformly distributed between [0,p] where p < 1 and the 
ones for the z-direction are in [0,1]. Next, a minimum 
random number is selected among the 2L random num- 
bers of the entire system, through which we determine 
the bead and the direction to move. Then the bead hav- 
ing the minimum random number is updated by shifting 
its location by the lattice constant a along the direction 
chosen. Next, the avalanche process of updating may oc- 
cur on neighboring beads when the separation between 
the nearest neighboring beads along the string becomes 
larger than In that case, the nearest neighbor- 

ing bead is also shifted by a lattice constant along the 
direction already selected through the minimum random 
number. The avalanche rule is then applied successively 
to next neighboring beads to hold the conservation of 
the separation between nearest neighboring beads. The 
dynamic rule we used is similar to the Sneppen dynamic 
rule but that the updating can occur in two different 
directions. Accordingly, we may say that our model is a 



vector Sneppen model. 
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FIG. 1. The initial flat configuration of the discrete ver- 
sion of elastic string. Each bead has two-component noises 
representing random pinning forces in y and z directions. 

The motion of the elastic string might be written by 
the following continuum equations, 

dty{x,t) = i^ydly + ^{d.^yf + r]y{x,y, z), 

dtz{x,t) = i^zdlz + y(3;j.z)^ +'qz{x,y,z), (5) 

which are coupled via the quenched noise. The justifica- 
tion of the above equations associated with the stochastic 
model is based on that our model is similar to the Snep- 
pen model but having the two different directions and 
the Sneppen model is a self-organized stochastic model 
belonging to the QKPZ universality in 1-1-1 dimensions. 
Even though the Sneppen model does not exhibit the PD 
transition, but the model gives the correct values of the 
roughness and the dynamic exponent in the QKPZ uni- 
versality. Accordingly, we think that the self-organized 
stochastic model we introduced in this paper could give 
the correct values of the roughness and the dynamic ex- 
ponents for the coupled QKPZ equations, Eq.(5). 

When p < 1, the beads are more likely to advance in 
the y-direction, because random numbers are accumu- 
lated in the smaller interval [0,p]. The ratio of the aver- 
aged velocities Vz/vy of the z and y directions is equal to 
p. Thus we can define the angle i/' — tan~^ p as the angle 
between the y-axis and the averaged propagating direc- 
tion. Transforming the coordinate system {y, z) by the 
angle ip into the new system (r|| ,r±), the above equations 
can be rewritten in terms of r|| and r±, where r|| is along 
the direction of the averaged velocity and r± is perpen- 
dicular to ry. The continuum equations in the (r\\,r±) 
coordinate system are written as 

dtr\\ix,t)= i^ndl^+i^i2dlr^ + ^{d^^Y + 



2 



dtrA^{x,t)^ V2idlr\\ + U22dlr^ + -^{d^r^f + (6) 
^{d^r\(f + ^((9x''x)(<9a;r||) + 77x(x,r||,r^), 

where Vab with a, 6 e 1 or 2 are functions of Vy, Vz, cosip, 
and sinV', and A^b with a,b G 1, 2, or 3 are of Ay, Az, cos ip 
and sin'0. 77 1| and ri± are the random noises in the longi- 
tudinal and the transversal direction, respectively. When 
p — 1, the advance of the elastic string in y and z direc- 
tions are identical, and the above equations are reduced 
simply as 

dt^{x,t) = v,,dl^ + ^(a.r||)2 -t- ^{d.r^f + r,||, (7) 
dtr^{x,t) = V22dlrx_ + ^{dxri_){dxr\0 +-q^. (8) 

Note that An 7^ in Eq. (7) and Eq.(8) is invariant 
under r±_ — In order to obtain the roughness ex- 
ponents for each direction, we consider the spatial corre- 
lation functions C\\ and C±_ after saturation, 

xi 

C^{x,t)^{^Y. (''^ i^ + ^i^t)- r± {^1 , t))Y', (9) 

xi 

which behave as C||(a:) ~ and C±{x) ^ x"^. Next, 
in order to obtain the growth exponents, we consider the 
temporal correlation functions C|| and Cj_, 

C\\{t2 - h) = {j ^2) - nit2) ~ r||(x,ii) + 

X 

C±{t2-h) = {jY.^r^{x,t2) - fAh) - r^{x,ti) + 

X 

r±{h)ry/', (10) 

where ti is taken as a saturation time. The correlation 
functions behave as Cy ^ (i2— ti)''" and C'_l ^ {t2—ti)^^ 
Numerical simulations are performed for the cases of 
p = 1 and p — 1/2. 

For p = 1, where Eqs.(7-8) are hold, the roughness 
exponents are measured as a|| « 0.63 and a± « 0.50, 
and the data are shown in Fig. 2. The growth expo- 
nents are measured as /3|| « 0.64 and i3± ~ 0.53, and 
the data are shown in Fig. 3. From the numerical results, 
the dynamic exponents are obtained as Z|| w 0.99 and 
z± ~ 0.95, which suggest that the true dynamic expo- 
nents be Zj| = zj^ = 1. This result is attributed to that 
the coherent effect propagates along the string and the 
chemical distance between two points on string remains 
invariant through the process of the dynamics . 
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FIG. 2. Plot of C|| (Cx) versus x for p = 1 in double loga- 
rithmic scales. The simulations are performed for the system 
size L — 2048. The lines obtained from the least square fits 
have the slopes 1.26 (top) and 1.00 
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FIG. 3. Plot of C|| {C±) versus time for p = 1 in double 
logarithmic scales. The simulations are performed for the sys- 
tem size L = 2048. The lines obtained from the least square 
fit have the slopes 1.27 (top) and 1.05 

Next for p = 1/2, the two spatial correlation func- 
tions for each direction are less distinctive than the case 
of p = 1 as shown in Fig. 4. The roughness exponents 
are obtained as a\\ « 0.60 and a± w 0.54. For the 
growth exponents, the time correlation function C|| in 
the longitudinal direction exhibits the power-law depen- 
dent behavior against time with the exponent /3|| « 0.64. 
However, for the transversal direction, the data do not 
exhibit a simple power-law type behavior, rather they 
show a crossover behavior from P± « 0.50 to (3± < 0.50 
in Fig. 5. From the numerical results, it is likely that the 
roughness exponents are ay = 0.63 and a± = 0.5, and the 
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dynamic exponents are z\\ = = \ even for < p < 1. 
Note that when p = 0, the dynamics in the transversal 
direction does not exist at aU in our modcL Accord- 
ingly, the numerical results indicate that the roughness 
and dynamic exponents do not depend on the angle ip. 
Also the numerical values of the roughness and dynamic 
exponents suggest that the dynamics in the longitudinal 
direction belongs to the directed percolation depinning 
universality class (the QKPZ universality class) in 1+1 
directions. On the other hand, for the transversal direc- 
tion, the numerical values a±_ w 0.5 and 2;x ~ 1 coincide 
with the ones in the anisotropic QKPZ universality class 
[10,14], even though the dynamics equation, Eq.(8), does 
not include the symmetric breaking term, (dxr±f'. We 
think that this coincidence is due to the presence of the 
term proportional to S^rj, in Eq.(8). 
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FIG. 4. Plot of C|| (Cx) versus a; for p = 1/2 in double log- 
arithmic scales. The simulations are performed for the system 
size L = 2048. The lines obtained from the legist square fit 
have the slopes 1.20 (top) and 1.08. 








FIG. 5. Plot of C\\ {C±) versus time for p = 1/2 in dou- 
ble logarithmic scales. The simulations are performed for the 
system size L = 2048. The lines have the slopes 1.28 (top) 
and 1.00 

In siimmary, we have introduced a stochastic model 
associated with the dynamics of the flux line in quenched 
media at the depinning threshold. The conservation of 
the total length of string through the dynamics in the 
current stochastic model leads to z = 1 in both the lon- 
gitudinal and the transversal direction. The roughness 
exponents were obtained as a\\ « 0.63 and a±_ ~ 0.5, 
respectively. 
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